AN INTEGRAL FORM OF THE NONLINEAR SCHRODINGER EQUATION 

WITH VARIABLE COEFFICIENTS 
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Abstract. We discuss an integral form of the Cauchy initial value problem for the nonlinear 
Schrodinger equation with variable coefficients. Some special and limiting cases are outlined. 



1. Introduction 



In the previous Letter we have constructed the time evolution operator for the linear one- 
Q_i. dimensional time-dependent Schrodinger equation 

zH^=H(t)^ (1.1) 

in a general case when the Hamiltonian is an arbitrary quadratic form of the operator of coordinate 
and the operator of linear momentum; see |15j . In this approach, several exactly solvable models, 
that have been studied elsewhere, are classified in terms of elementary solutions of certain charac- 
^ ■ teristic equation. A particular solution of the corresponding nonlinear Schrodinger equation with 
variable coefficients had been obtained in a similar fashion. In this paper we rewrite a nonlinear 
g ! Schrodinger equation 

g : (*§i~ H ^ ^ (x > t) = F (t ' x ' ^ (x ' t)] (1 - 2) 

o 



in an integral form and consider several examples. In general, we do not assume that the time- 
dependent linear Hamiltonian H (t) here is the Hermitian operator. 

2. A General Lemma: Duhamel's Principle 

The following result is helpful in study of solutions of the nonlinear Schrodinger equation (11.21) 
by a fixed point argument. 

Lemma 1. Suppose that the Cauchy initial value problem 

l Wt~ H (t) ) ^° (x ' f) = °' ^° (x ' t} '*=° = v (x) (2,1) 

for a linear time- dependent Schrodinger equation can be solved in terms of the time evolution oper- 
ator 

rl> (x,t) = U(t)if> (x,0) (2.2) 
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with 

(i^-H(f)W) = 

and 

U(0) = I = U(t)U(ty 1 . 

Then the initial value problem 

*j t -H (t)j if) = F(t, x, if)) , if) (x, t) | t=0 = X (x) 

for the nonlinear Schrddinger equation can be rewritten as an integral equation 

ifj (x, t) = U (t) if) (x, 0) - % [ U (t)U (s)' 1 F(s,x,if)(x,s)) ds 

Jo 

in terms of the time evolution operator U (t) for the linear equation and its inverse. 
Proof. Indeed, from (12. 6p 

- H (t)j if)=U^-H (t)j U (t) if) (x, 0) 



+i / H (t)U (t)U (s)" 1 F (s,x,ip (x,s)) ds, 



where 



^ J C/(0C/(s) _1 F( S ,x^(x,s)) ds 
= U(t)U(ty 1 F(t,x,if>(x,t)) 

^lu{sy l F{s,x^{x,s)) ds. 



Thus, by Q and fl2T4j) 



i^-H(t)jif) = F(t,x,if>(x,t)) 

-ij ^ t -H{t)^jU{t){U(sy l F(s,x,iP(x,s))) ds 
= F{t,x,if){x,t)). 
Initial conditions are satisfied in view of 

lim f UtyUisy 1 F(s,x,ip(x,s)) ds 



U(0) lim / Uis)' 1 F{s,x,ip{x,s)) ds = 0, 
*^ 0+ Jo 
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where we have used the fact that if lim^ + /(£) = / (0) , then 

lim I f(s) ds = 0. 
t->o+ Jo 

This completes the proof. □ 

When F does not depend on ip, expression (12.60 solves the Cauchy initial value problem (12.51) for 
the corresponding nonhomogeneous Schrodinger equation. 

The integral equation (12.61) can also be rewritten in the form 

V>(x,0) = U- 1 (t)^(x,t)+t [ Uis)' 1 F{s,x,tfj{x,s)) ds (2.7) 

Jo 

= T(t)- X ij)(x,t), 

which gives explicitly the inverse of the time evolution operator 

ip (x, t) =T (t) ip (x, 0) (2.8) 

for the nonlinear Schrodinger equation (11.21) . Thus, in general, solution of the Cauchy initial value 
problem ( 12.51) is a problem of inversion of this nonlinear integral operator. 

3. Quadratic Hamiltonians 

The fundamental solution of the linear Schrodinger equation with the quadratic Hamiltonian of 
the form 

^ = -a (*) ^ + b (t) *V - i [c (t) x^ +d(t)^-f (f) xij + ig (t) (3.1) 

where a(t) , b (i) , c(t) , d(t) , / (t) , and g (t) are given real- valued functions of time t only, can be 
found with the help of a familiar substitution 

$ = Ae lS = A (t) e iS{x ^\ A = A (*) = 1 (3.2) 

y2irin (t) 

with 

S = S (x,y,t) = a (t) x 2 + (5 (t) xy + 7 (t) y 2 + 5 (t) x + e (t) y + k (t) , (3.3) 

where at(t) , (5 (t) , 7 (t) , 8 (t) , e (t) , and k (t) are differentiable real- valued functions of time t only; 
see [15]. Indeed, 

dS fdS\ 2 , 2 , , ,dS 



at -"\T X ) -fa J + /x + fa-«) s (3.4) 

provided 

— = a—- + d = 2a(t)a(t) + d(t) . (3.5) 
2fi ox 2 

Equating the coefficients of all admissible powers of x m y n with < m + n < 2, gives the following 
system of ordinary differential equations 

dot 

— + b(t) + 2c(t)a + 4a(t)a 2 = 0, (3.6) 
dt 

^ + {c{t)+4a{t)a{t))/3 = 0, (3.7) 
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^ + a(t)(3 2 (t) = 0, (3.8) 
dS 

— + (c (t) + Aa (t) a (t)) 5 = fit) + 2a (t) g (t) , (3.9) 

de 

— = (g(t)-2a(t)6(t))P(t), (3.10) 

^- = g (t)6(t)-a(t)6 2 (t), (3.11) 

where the first equation is the familiar Riccati nonlinear differential equation; see, for example, |23j, 
[5%] . [4"5] and references therein. Substitution of (I3.5P into (I3.6P results in the second order linear 
equation 

n"-r(t)n' + 4ff(t)ii = Q (3.12) 
with ^ f 

T (t) = — - 2c + 4d, a{t)=ab-cd + d 2 + -(- -^-\ , (3.13) 
a 2 \ a a J 

which must be solved subject to the initial data 

/i(0) = 0, n'(0) = 2a (0)^0 (3.14) 

in order to satisfy the initial condition for the corresponding Green function; see the asymptotic 
formula (13.221) below. We refer to equation (I3.12p as the characteristic equation and its solution 
f/,(t), subject to (I3.14p . as the characteristic function. As the special case (13.121) contains the 
generalized equation of hypergeometric type, whose solutions are studied in detail in [35]; see also 
®, m, m, and @5]. 

Thus, the Green function (fundamental solution or propagator) is given in terms of the charac- 
teristic function 

ip = G (x, y, t) = 1 e i(a(t)x 2 +{3(t)xy+7(t)y 2 +mx+e(t)y+<t)) _ (3.15) 
y^2nifi (t) 



Here 



a(t)= j-m_m (3 . 16 ) 

V ; 4a (t) n (t) 2a (t) ' V ; 



1 < '* 



7(*) = T^L exp( -2 / (r(7-)-2r/(7)) r/7-1 (3.18) 



w (t)/l' (t) 

- 4 / , ^^(- p (- 2 r (c(A) - Mw),tt )) ,fr " 

5(t) = -^—exp(-[ (c(T)-2d(r)) dr) (3.19) 



x / exp ( / (c(A) - 2d (A)) dA 

^(r))/»(r)H _ , , 
a (r) / 2a (r) 



x ( ( / (r) - ^4 (r)") /, (r) + (r) ) dr, 
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£« = ~jM S (t) expf- j\c(r)-2d(T)) dr\ (3.20) 
+8jf ^^exp(-jf (c(A)-2«i(A)) d\j (/i(t)*(t)) <ir 

+2 jf ^7 exp (" I (e (A) - 2<i (A)) dX ) ( f (r) - ^ (T) ) dT ' 

* M = ^ (*) - 4 /' 0. W * M)' * (3.21) 
V'[t) Jo (jjl'(t)) 

' " {T) ^(r)6(r))(f(r)- d -^-g(r)) dr 



/ //(r) V a ( T )' 

with 5 (0) = g (0) / (2a (0)) , e (0) = — 5 (0) , and ft (0) = 0. Integration by parts has been used to 
resolve the singularities of the initial data. Then the corresponding asymptotic formula is 

G (x, y, t) = -gj ~ 7 =L= exp (i^gj exp (, - „)) (3.22, 

as t — > + . Notice that the first term on the right hand side is a familiar free particle propagator 
(cf. (E3D below). 

By the superposition principle, an explicit solution of the Cauchy initial value problem 

i^ = H(t)1>, i>(x,t)\ t=0 = p(x) (3.23) 

on the infinite interval — oo < x < oo with the general quadratic Hamiltonian as in (13.11) has the 
form 

/oo 
G{x,y,t) V(2/,0) dy. (3.24) 
-oo 

This yields the time evolution operator ([2.2]) explicitly as an integral operator. 

4. Inverses of the Time Evolution Operator 

In the previous section we have discussed how to construct the time evolution operator for the 
linear Schrodinger equation with the quadratic Hamiltonian (13.11) . In this section we study the 
inverses 

/oo 
G(x,y,t) <%,0) dy, (4.1) 
-oo 

/oo 
H(x,y,t) ^(y,t) dy (4.2) 
-oo 

such that 

U it) U- 1 (t) = U~ l (t) U(t) = I = id. (4.3) 

Here we introduce 

p iS(x,y,t) 
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H(x,y,t) = G*(y,x,t)esxp(- J (c(r)-2d(r)) dr^ (4.5) 

-i,S(y,x,t) / ft \ 

= 7^m^\~L (c{T) - 2i{T))dT ) 

with S (x, y,t) = a (t) x 2 + [3 (t) xy + >y (t) y 2 + 8 (t) x + e (t) y + k (t) . 

First we formally prove the following orthogonality relations of the kernels 

G(x,y,t)H(y,z,t) dy = e *H*) § ( x _ z ) ? ( 4 . 6 ) 
H(x,y,t)G(y,z,t) dy = e ^(7(t)(*+*)+*(t))(*-*) § ( x _ z ) ? (4.7) 



-00 
00 



where 5 (x) is the Dirac delta function with respect to the space coordinates (do not confuse with 
a given function of time 8 (t) throughout the paper). 



1 ' rt 



Indeed, by (|Pj) - (|3~5j) one gets 

G(x,y,t)H(y,z,t) dy = — ^— exp ( - / (c(r)-2d(r)) dr 
> M f J \ Jo 

xe i(a(t)(x+z)+6(t))(x-z) J_ / e i(3{t){x-z)y dy 
= e i(a(t)(x+z)+5(t))(a;-z) £ ^ _ ^ 

in view of (I3.17P with —(3 (t) > and the integral 



00 



H0 = ir I ^ dt (4. 



2tt 

as the Dirac delta function. The formal proof of ( 14. 7ft is similar and is left to the reader. 
Now we have 

CT 1 (t) C/ {t) V> (x, 0) = LT 1 (t) V> (x, t) 
H(x,y,t) ip(y,t) dy 



■00 
00 



■00 

00 / /"OO 



00 \J —00 

00 



H(x,y,t) (J G(y,z,t) iP(z,0) dz j dy 
H (x,y,t)G(y,z,t) dy^j ip{z,0) dz 

iW)(*+*)+e<t))(*-*) S(x-Z) if) (*, 0) dz 

= VM), 

or (t) [/ (t) = I. A formal proof of the second relation U (t) U^ 1 (t) = I is similar and left to 
the reader. 
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An integral form (12. 6p of the nonlinear Schrodinger equation (jl.2p contains also the integral 
operator U (t, s) = U (t) U' 1 (s) : 



U (t) U- 1 (s) if>(x,s)= I G (x, y, t, s) (y, s) dy (4.9) 

with the kernel given by 

G(x,y,t,s)= I G(x,z,t)H(z,y,s) dz. (4.10) 



Here 

f°° 1 

/ G (x,z,t) H (z,y,s) dz = -^=^= (4.11) 

J-oo 2n ^ n {t) n (s) 

xexp^-^ (c(r)-2d(r)) dr^j 

xe i(a(t)x 2 -a(s)y 2 +8(t)x-8(s)y+K(t)-K(s)) 

X I™ ^{(^-l^ + m^-my+e^-eis^z) dz 
J — oo 

and with the help of the familiar elementary integral 

J(az*+2bz) dz = /Z^ e -*> 2 /» (4 . 12 ) 



oo 

f 

V a 



see Refs. [3] and [36], we get 



G (x ' *> *' s) = i a ■ m ? w n ^ (4 - 13) 

y/4mfi (t)//(s) (7(5) -7^)) 



x exp ^— y (c (r) - 2d (r)) dr 

x exp (i (a (t) x 2 — a (s) y 2 + 5 (t) x — 5 (s) y + n (t) — k (s))) 
'(/?(t)x-/?( S )y + £ (t)-£( S )) 5 



x exp 



4i ( 7 (t)-7(s)) 



This can be transform into a somewhat more convenient form 

G(x,y,t,s)= , 1 expf-/ (c (r) - 2d (r)) dr ] (4.14) 



x exp 
x exp 
x exp 



( g (t)- £ ( g )) 2 -4(7(t)-7( g ))(^(t)-^( g ) ) 
4z (7 (t)- 7 (s)) 

(e (t) - e (g)) (g (fls-j (g) y) - 2 ( 7 (t) - 7 (fj) (j (t) x - j (g) y) 

2, ( 7 (t)- T (s)) 

(g (t) x - g (g) y) 2 - 4 (7 (t) - 7 (a)) (a (t) x 2 - a (a) y 2 ) ' 
4, ( 7 (t)- 7 ( s )) 



In the limit s — > t with s < t one arrives at the kernel of the identity operator. We leave the details 
to the reader. 
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5. Axillary Tools: An Estimate and a Functional Equation 

Consider the operator U (t,s) = U (t) U^ 1 (s) . From (14.91) and ( 14. 14j) one gets 

\U(t,s)4(x,s)\ = ! 

J — oo 
oo 



< 



G(x,y,t,s)ip(y,s) dy 

\G(x,y,t,s)ip(y,s)\ dy 
1 

V47T|/x(t)/x( S ) ( 7 ( S )- 7 (t))| 

x exp ^— J (c(r) — 2d(r)) dr 

POO 

x / \*l>{y,s)\ dy 



and as a result 



\U{t,s)il>(x,s)\ < 



y/4Tr\ti(t)fj,(s) (7(s)-7(*))l 



x exp 



(c(r) -2d(r)) dr 



x / hMy s s)l 



Thus the familiar estimate 
\\U{t,s) 



< 



y/4n\fj,(t)ii(s) (7 (s)-l(t))\ 



x exp 



(c(r)-2d(r)) dr 



(5.1) 



(5.2) 



holds in the case of the general quadratic Hamiltonian (13.11) (cf. [26J). 

As we shall see in the next section, some solutions of the characteristic equation (13.121) obey the 
following property 

( t + s 

fi(t)(i(s)(j(a)-'y(t))=x 



Then 



„( ()/J ( s )!M^W = , 



t - s 



jj, (t — s) . 



t + s\ n (t — s) — // (0) 
t-s 



(5.3) 



and in the limit s — ► t one gets 
or by (13T4D 



c?7 



V (*)^ = x(*V (0) 



dt +2a(0) ^(r) =0 - 



(5.4) 
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But, in view of $£ffy and (j3T7|) . 



dl + exp [ -2 / (c (r) - Id (r)) rfr I -^TTV = °- ( 5 - 5 ) 



\ o / 
Therefore, the addition property (15.31) may hold only when 

X (t) = ~ exp ^-2 j (c (r) - 2d (r)) dr J (5.6) 

Some examples will be given in the next section. 

6. Examples 

Now let us consider several elementary solutions of the characteristic equation (13.1 2ft ; more com- 
plicated cases may include special functions, like Bessel, hypergeometric or elliptic functions [lj, 
[55] . [37) . and [IS]. Among special cases of general expressions for the Green function fl3.15l) - (l3.21l) 
are the following 



6.1. A Free Particle. When a = 1/2, b = c = d = f = g = 0, and jj" = 0, [i = t, one gets 

^•'^"Kt) (6 ' 1) 

as the free particle propagator [21]. In this case a = —(3/2 = 7 = 1/ (2t) and an elementary identity 

(x/t-y/sf x 2 y 2 _ (x-yf 



(1/t-l/s) t s t-s 
implies that 

G (x, y, t, s) = G (x, y, t - s) (6.2) 

from the general formula (I4.14p . The time evolution operator of the linear problem is given explicitly 
as the following integral operator 

U (t) X (x) = -== / e^~y?/ 2t x (V) dy (6.3) 

\J2l\lt J~oo 

and the traditional nonlinear Schrodinger equation 

d 1 d 2 \ 

*V + 1 ^ = AIV'I 2 ^ ib, A = constant, < v < 1 (6.4) 

at 2 ox J 

has the familiar integral form 

$ (x, t) = U (t) ib (x, 0) - iX [ U{t-s) \ip {x, s) | 2 *> (x, s) ds. (6.5) 

Jo 

See [14] . [16] . [26] . [H] and references therein for more information. 
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6.2. A Particle in a Uniform External Field. For a particle in a constant external field, where 
a = 1/2, b = c = d = g = and / = constant, fi = t, the propagator of the linear problem is given 
by 

i(x-yf\ fif(x + y) if 2 



G {x - y - t) = 7m exp ) exp {^r^t - i?) ■ (e 6) 

This case was studied in detail in [2J, [B], [2TJ, [21], [33] and [40]. We have corrected a typo in |21j; 
see [H] for a complete list of known errata in the Feynman and Hibbs book. 

In this case once again a = —(3/2 = 7 = 1/ (2t) and, in addition to the case of a free particle, 
5 = s = (ft) i '2 and k = —f 2 t 3 /24. An elementary calculation shows from the general expression 
(14.141) that relation (16.21) holds. Therefore, the corresponding nonlinear Schrodinger equation 

d Id 2 „ \ . 



i— + 2~q~^2 ~^ f x J ^ = ^ 1^1 V ^' ^ = cons t an t, < i/ < 1 (6.7) 

has the integral form (16.51) . where the linear propagator is given by (16.61) . 

In a more general case of a particle in a uniform electric field changing in time with a similar 
velocity-dependent term 

d . i d 2 . „ a \ , . , , |2v 



i to + 2W +nt)x - i9 ®te** = m (6 ' 8) 



where / (i) and g (t) are functions of time only, the propagator of the linear problem has the form 

1 I i(x — 



G(x,y,t) = -jL= exp ( ^ — ^- I exp(i($(t)ar + e(t)2/ + «(t))) (6.9) 



with 



= 7 / (/ (r) r + g (r)) dr, (6.10) 
e(t) = -<?(*) + / f(r) dr (6.11) 



and 



t 



K(t) = -5 2 (t)- I r5(t)f(r) dr. (6.12) 



A semigroup property [T3], related to (16.21) . does not hold anymore and one has to use a general 
expression (14.141) in order to write the integral equation (I2.6p . But, in view of an elementary identity 

H (t) fi (s) ( 7 (s) - 7 (*)) = ±(Ji(t-s), (6.13) 

an important addition formula still holds for the amplitude of the kernel G (x, y, t, s) in (I4.9I) - (I4.10I) 
and fl4TT4"j) of the operator U (t, s) = U (t) U~ l (s) . 
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6.3. The Forced Harmonic Oscillator. The simple harmonic oscillator with a = b = 1/2, 
c — d — f — g — and fi" + fi = 0, fi = sint has the familiar propagator of the form 

G (x, y, t) = 1 exp ( -4— ( (x 2 + y 2 ) cos t - 2xy) ) , (6.14) 

which is studied in detail at @], [22], [25], [30], [32], 03]. 

Once again relation (16. 2p holds and the corresponding nonlinear Schrodinger equation 
dib ltd 2 \ 

i-zr + -hr^-£ 2 W' = A U\ 2v ib, A = constant, < v < 1 (6.15) 
at 2 / 

has the integral form (16.51) . where the propagator is given by (16.141) . See [H], [8], [9], [TO], [H], [12] . 
[T3] and references therein for investigation of solutions of this integral equation by a fixed point 
argument. 

A linear problem extension to the case of the forced harmonic oscillator including an extra 
velocity-dependent term and a time-dependent frequency is discussed in [IT], [18], [21] and [29J. 
The nonlinear Schrodinger equation of interest is 

i^ + l -aAij> + f (t) xi> -ig(t)^ = \ \ih\ 2u ij) (6.16) 



dt 2 \dx 2 J dx 
and the corresponding propagator has the form 

G(x,y,t) = 1 . exp(-^—((x 2 + y 2 )cost-2xy)) (6.17) 
x exp (i (S (t)x + e(t)y + k (t))) 



with 



1 f l 

5{t) = - — / (/(r)sinr + g(r)cosr) dr, (6.18) 



sint J 



6(t) rsmrS(T) , rfir) , 

£ (t) = 4+ / ^dr + / dr 6.19 

cost ./o cos z r J cost 



and 



«(*) 



1 1 f 

- tant 5 2 (t)-- tan 2 r <5 2 (r) rfr (6.20) 

2 2 Jo 

— / tanr 5 (t) f (t) dr. 
Jo 

The addition property (15.31) holds in this case. We leave the detail to the reader. 



6.4. A Modified Oscillator. Furthermore, an exact solution of the n-dimensional time-dependent 
Schrodinger equation for certain modified oscillator is found in [31] . In the one- dimensional case we 
get functions a = cos 2 1, b = sin 2 1, c = 2d = sin 2t and our characteristic equation (13.121) takes the 
form 

// + 2tant //-2// = 0, (6.21) 
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whose elementary solution \i = cos t sinh t + sin t cosh t satisfies the initial conditions (13.141) . Further, 
the corresponding propagator is given by 

G (x, y, t) = 1 = (6.22) 
a/ 2ni (cos t sinh t + sin t cosh t) 



x exp 



; x 



2 „,2 



y 2 ) sin t sinh t + 2xy — {x 2 + y 2 ) cos t cosh t 



2i (cos t sinh t + sin t cosh t) 

which was found in [31] as the special case n — 1 of a general n-dimensional expansion of the 
Green function in hyperspherical harmonics. We have showed that f)6.22p is a generalization of the 
propagator for the simple harmonic oscillator; see Ref. [31] for more details. 

The corresponding nonlinear Schrodinger equation [15] 



M , *t&ll> ■„■>, ,1,, , ; , A,./ 7 ''' , ,,\ _,./,M.,l->" 



i— — h cos t — — r- — sin t x ip + i sin t cost ( 2x— — \-ip ) = h (t) \ip\ v ip (6.23) 
at ox 2 \ ox J 

can be rewritten in an integral form. We leave the detail to the reader. 
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